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Retrieval of classical behaviour in quantum cosmology is usually discussed in the framework 
of rmcfeuperspace models in the presence of scalar fields and the inhomogeneous modes corre- 
sponding either to gravitational or scalar fields. In this work, we propose an alternative model 
to study the decoherence of homogeneous and isotropic geometries where the scalar field is 
replaced by a massive vector field with a global internal symmetry. We study here the cases 
with U(l) and 50(3) global internal symmetries. The presence of a mass term breaks the 
conformal invariance and allows for the longitudinal modes of the spin-1 field to be present 
in the Wheeler-DeWitt equation. In the case of the U(l) global internal symmetry, we have 
only one single "classical" degree of freedom while in the case of the SO (3) global symmetry, 
we are led to consider a simple two-dimensional minisuperspace model. These minisuperspaces 
are shown to be equivalent to a set of coupled harmonic oscillators where the kinetic term of 
the longitudinal modes has a coefficient proportional to the inverse of the scale factor. The 
conditions for a suitable decoherence process and correlations between coordinate and momenta 
are established. The validity of the semi-classical Einstein equations when massive vector fields 



(Abelian and non-Abelian) are present is also discussed. 
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1 Introduction 



The emergence of classical properties from the quantum mechanics formalism is still largely 
an open problem. Some progress has, however, been achieved through the so-called decoherence 
approach. On fairly general grounds, the decoherence process takes place as one considers the 
system under study to be part of a more complex world that interacts with other subsystems, 
usually referred to as the "environment" . The latter usually consists in the set of unobserved 
or irrelevant degrees of freedom which are traced out, implying that, at least in an operational 
way, the wave-function evolves non-linearly and is lead to its collapse |TJ . In this way, quantum 
interference effects among states of the system are suppressed by the interaction with the 
environment. This suppression comes about as one integrates out the irrelevant degrees of 
freedom. This coarse-graining procedure leads to an effective action of the original model and 
clearly generalizes the usual field fluctuation splitting that is adopted when carrying out a 
quantum loop expansion in the standard effective action calculation. From the operational 
point of view, decoherence in the context of quantum cosmology, implies that the Universe is 
essentially an open system as observers have necessarily to disregard large classes of variables 
in any relevant observation. In addition to the notion of decoherence, which strictly means that 
there is no interference among different quantum states, a further condition that a system must 
satisfy to be regarded as classical is, of course, that it is driven by classical laws, implying that 
a sharp correlation between configuration space coordinates and conjugate momenta should 
exist in the wave function. 

These ideas have been developed in some depth in the context of quantum cosmological 
models where considering the decoherence process is mandatory as the observable 



Universe behaves clearly in a classical fashion and one expects that this classical features 
arise from a more fundamental quantum description. Furthermore, although minisuperspace 
models can be justified on symmetry grounds, the truncation which turns the full quantum 
gravity problem with its infinite degrees of freedom into a problem with a finite number of 
degrees of freedom actually violates the uncertainty principle, as the amplitudes and momenta 
of inhomogeneous modes are set to zero and the non-linear interactions of those modes among 
themselves and with the minisuperspace degrees of freedom are disregarded. Moreover, the 
validity of the minisuperspace approximation (cf. || and references therein) is ensured only 
when the back-reaction of the inhomogeneous modes on the minisuperspace variables is shown 
to be small. 

Most of the literature concerning the emergence of classical behaviour in the context of 
quantum cosmology considers scalar fields Pfl-flOfl, where the environment corresponds to the 



inhomogeneous modes either of the gravitational or of the scalar fields M- \ 1 1]. This is a short- 



coming of most of the decohering models discussed so far as prior to the inflationary epoch the 
Universe was dominated by radiation, i.e. Yang-Mills fields, and, after spontaneous symmetry 
breaking phenomena, by massive vector fields. Furthermore, a relevant issue concerning the 
decoherence approach is to stablish to what extent some of its features are specifically related 
with scalar fields and to achieve this goal one has to consider other fields to play the role of 
environment. 

In this paper we propose a model to study the decoherence of homogeneous and isotropic 
geometries, where the scalar field is replaced by a massive vector field with a global internal 
symmetry. Our aim is to assess if such models are on the same footing with the ones where the 
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quantum to classical transition is achieved via tracing out the higher modes of scalar fields, or 
in other words, if massive vector fields are equally effective in playing the role of an environment 
with respect to the observed "classical" degrees of fredoom. Notice that the presence of a mass 
term is a crucial feature as it breaks the conformal symmetry of the spin-1 field action which 
leads to a Wheeler-DeWitt equation where the gravitational and matter degrees of freedom 
decouple [12], similarly to the case of a free massless conformally invariant scalar field []13" 



Hence, the presence of a mass term allows for interaction between gravitational and matter 
degrees of freedom providing a scale at which the decoherence process can take place. Moreover, 
the breaking of conformal invariance makes it possible for the longitudinal modes of the spin-1 
field to be present in the Wheeler-DeWitt equation. As we shall see, the resulting system 
is equivalent to a set of coupled harmonic oscillators where the kinetic energy terms of the 
longitudinal modes have a coefficient proportional to the inverse of the scale factor. 

An important point related to the discussion of quantum cosmology with vector fields con- 
cerns the compatibility of the simple homogeneous and isotropic Friedmann- Robert son- Walker 
(FRW) geometries we shall consider and the matter content of the Universe. As far as the 
homogeneous modes of the vector fields are concerned, we shall use here the Ansatz formulated 
in Refs. |L^], |pl}]-[P0|] . More specifically, since the physical observables have to be SO (4)- 
invariant, the isometry group of closed FRW cosmologies, the fields with internal degrees of 
freedom may be allowed to transform under SO (4) if these transformations can be compensated 
by an internal symmetry transformation. Fortunately, there is a large class of fields satisfying 
the above conditions, namely the so-called SO(4) -symmetric fields, i.e. fields which are invari- 
ant up to an internal symmetry transformation. However, such construction can only be made 
possible in the presence of non-Abelian symmetries and hence, for the the Abelian case, we 
shall require the homogeneous modes of the spatial components of the vector field to vanish. 
For a R x S 3 topology, such minisuperspace constructions with homogeneous and isotropic 
metric and gauge fields with SO(N) and SU(M) (N > 2, M > 3) gauge groups have been 
recently carried out. For Euclidean FRW geometries Einstein- Yang- Mills wormhole solutions 
have been found in Ref. |]I4" |. Classical Einstein- Yang- Mills solutions for closed geometries 
and SO(N) gauge groups were obtained in Ref. |15|| . The ground-state wave-function of the 
Einstein- Yang-Mills system with SO(3) gauge group was found in Ref. [|12|]. Massive vector 



fields with SO (3) global symmetry in flat FRW universes were studied in Ref. [16]. Finally, 
in a R 4 topology, homogeneous and isotropic metric and Yang-Mills field configurations were 



considered in Ref. [17] and the inclusion of the dilaton in the context of string theories was 



carried out in Ref. 18 



In this work, we shall study the decoherence process for homogeneous and isotropic metrics 
in the presence of massive vector fields with U(l) and 5*0(3) global symmetries in a R x S 3 
topology. For simplicity, we begin by considering the case with U(l) global symmetry. The 
role of the environment will be played by the inhomogeneous modes of the U(l) field and the 
corresponding minisuperspace will be actually one-dimensional as it has only a single "classical" 
degree of freedom as physical observable, the scale factor. Although this choice of the matter 
content may seem rather restrictive, we shall see that some of our results may be extended to the 
richer and more interesting system with S0(3) global symmetry, which we analyse afterwards. 
In particular, we shall establish, for both cases, the conditions required for achieving correlation 
between coordinates and momenta and a satisfactory decoherence process. The validity of the 
semi-classical Einstein equations when massive vector fields are present will be also discussed. 
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The minisuperspace for the non-Abelian model hereby studied is two-dimensional due to the 
specific Ansatz for the homogeneous modes of the spin-1 field (see Refs. ||14|-||20fl). 



Multi-dimensional minisuperspace models have of course, a much richer structure and are 
therefore far more interesting to consider in what regards the retrieval of classical behaviour 
[[7], ^], [M], In models containing one single "classical" degree of freedom, the Hamilton- 

Jacobi equation has only two solutions, generating the same trajectory in opposite directions. 
The semi-classical wave-function has two WKB components, each of which may be called a 
WKB branch, of the form 



9[0,£\ = J2C {n) [0}e iM ^°^ {n) [0,S}. (1) 

(n) 

Here the subindex (n) labels the WKB branches (taking only two values, say ±1, uniquely 
identifying the two possible solutions of the Hamilton- Jacobi equation) and O, £ denote, re- 
spectively, the "classical" physical observables and the extra degrees of freedom corresponding 
to the environment. To achieve a proper decoherence, one needs not only that the reduced den- 
sity matrix turns out to be diagonal, but also that the different WKB components in (|XJ) have 
negligible interference among the diagonal terms. It is important to stress that the analysis of 
correlations should be done within each classical WKB branch (i.e. a diagonal term, n = n'). 
The interference effects between the two possibilities of moving along the one- dimensional clas- 
sical trajectory (corresponding to the expanding and collapsing wave function components, 
respectively) were shown to be effectivelly suppressed ||], which is interpreted as particle cre- 
ation ||]. 

For a system with N degrees of freedom, the Hamilton- Jacobi equation is expected to have 
an N-l parameter family of solutions, each one generating a N-l parameter family of classical 
trajectories in the minisuperspace. In the multi-dimensional case, a general solution of the 
Wheeler-DeWitt equation may contain an infinite superposition of semi-classical solutions of the 
form ([I]) with the subindex (n) now corresponding to a set of parameters that uniquely identify 
a specific Hamilton- Jacobi solution. However, each WKB branch must actually be interpreted 
as describing a whole family of classical trajectories, i.e. a set of different universes and not a 
single one as for the N = 1 case. Furthermore, the N = 2 and N > 2 cases are rather different 
as far as the diagonalization of the reduced density matrix is concerned. These differences seem 
to have gone unnoticed until they were discussed in Ref. [Bf. An example of N = 2 model was 



studied in Ref. |24| , where the two "classical" degrees of freedom correspond to the scale factor 
and homogeneous mode of a minimally massless scalar field and the environment was identified 
with the inhomogeneous perturbations of another minimally massless scalar field. The multi- 
dimensional cases also allow one to better address the relation between the reduced density 
matrix formalism and the Feynman- Vernon influence functional E7J, and to the Schwinger- 



Keldish or closed time path effective action |28| , as pointed out in Ref. . 



Some preliminary work performed in Ref. [^] , with homogeneous and isotropic cosmologies 



in the presence of massive vector fields, has shown that the necessary ingredients for the process 
of decoherence to take place are present, although several aspects of the full hyperspherical 
harmonics expansion of the fields remained to be fully assessed. The Wheeler-DeWitt equation 



obtained in Ref. |L9[ is fairly similar to the one corresponding to a FRW minisuperspace 
model with a massive scalar field with a quartic self-coupling, A0 4 , conformally coupled with 
gravity M. In addition, we shall compare the results of this paper with the ones obtained 
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previously in Refs. [0, |8|]. Actually, there are many similarities and this will allows us to 
use some of the framework described in those references^. We would like to mention that, in 



considering perturbations in the quantum Einstein- Yang-Mills model of Ref. ||12j| , the authors 
of Ref. p2] consider also, as we do, the harmonic expansion of the fields involved on 5*0(4) 
and SO(4)/SO(3) (see above). 

This paper is organized as follows. In Section 2 we present our model and, through general 
considerations, introduce an Ansatz for the vector fields with global symmetry as well as the 
expansion in 5 3 hyperharmonics which will give rise to an effective action. We then proceed to 
a rmcfeuperspace description of such a model, which will allow us to study both the U(l) and 
50(3) global symmetry cases. In Section 3, we discuss the decoherence process and correlations 
within the U(l) model and, in Section 4, we address the perturbed minisuperspace model with 
massive vector fields with 50(3) global symmetry. Our conclusions are presented in Section 5. 



2 FRW Minisuperspace with Massive Vector Fields 



The action of our model consists of a Proca field coupled with gravity: 
5 = J d 4 x\^g 



±{R - 2A) + ±T T (F$F^) + 1 m 2 Tr(A^ A«") 



(2) 



where k 2 = SiiMp 2 , M P being the Planck mass, e is a gauge coupling constant and m the mass 
of the Proca field. To action (|^) one adds the boundary action 



S B = -— f d 3 xVhK 
k z JdM 



(3) 



with hij(i,j = 1,2,3) being the induced metric on the three-dimensional boundary dM of M, 
h = det(hij) and K = Kj^ is the trace of the second fundamental form on dM. 

In quantum cosmology one is concerned with spatially compact topologies and we will 
consider here the FRW Ansatz for the Rx S 3 geometry 



ds = a a (77) 



(4) 



wheref] a 2 = 2/3ir, rj is the conformal time, N(rj) and a(rj) being the lapse function and the 
scale factor, respectively and u % are left-invariant one-forms in 577(2) ~ 5 3 which satisfy 



duf 



-€kijUJ l A a/ . 



(5) 



Aiming to obtain solutions of the Wheeler-DeWitt equation satisfied by the wave function 



G. 



ijki 



S 2 M 2 

+ Mj > ^(^R-2A) + 1 ^VhT 



Sh^Sh 



5A { , a) 



nh t3 ,Af] = 0, (6) 



1 The authors are grateful to B.L. Hu and J. P. Paz for pointing this out. 

2 Usually one makes the choice tr 2 = 2/3irMp, but we shall keep the powers of Mp explicitely in our effective 
action in order to compare our results with the ones of Refs. M 0, yjl. 
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where the superspace metric is given by 



Gijki = (hikhjt + hithjk — hijhkg) , 



(7) 



we shall expand the metric as 



hij = a a (fiy + €ij) , (8) 

with fijj being the metric on the unit S 3 and a perturbation that can be expanded in scalar 
harmonics V J N M (g), which are the usual (2J+ l)-dimensional SU(2) matrix representation and 
spin-2 hyper spherical harmonics M N (g) on S 3 as: 



= ^ E \/^/m(^VG?) + <° n i (Vm n) *a 



(9) 



where 



e,4 



E 



32(n 2 - 4) 



\| 15(n 2 - 1 
+ E MdffWY, 

\J-L\=2 V 



MN 



(9)+ E 

|J-L]=1 



'32(rr 2 - 4) 



JMN ( 



2 LJ 
m MN 



(10) 



with cr™ described below, the coefficients aj , . . . , dj^ depend only on the conformal time, 
(n m m) represent 3-J symbols or Clebsh-Gordon coefficients and n = J + L + 1. 
The massive vector field 



Aab, .msr Aab^m, Arr 

A m uj s l ab = A m a i u J a b, 



111 



where denote the one- forms in a spherical basis with m = 0, ±1, T ab are the SO (3) group 
generators and 

" 



i 1 

— — 

i 



can be expanded in spin-1 hyperspherical harmonics as [f2~Tf : 



Ai(r],x J ) 



E 



a 



abJM 
N 



( V )V J N M (g)T ab = 0+ £ 



a 



JMiV J'M'N' 



abJ'M' 
N' 



'MV J ' N F{g)T abl (13) 



LJNM 
1 



, 2a 
3n 



(12) 



^ 6 T a6+ £ /S^fa) ^M^K™ %b, (14) 
L'J'N'M' 



where a = e 2 /47r and Ao is a scalar on each fixed time hypersurface, such that it can be expanded 
in scalar harmonics T? J N M (g). The coordinates x l are written as an element of 577(2). The 
expansion of A{ is performed in terms of the spin-1 spinor hyperspherical harmonics, mn{.9)- 
Longitudinal and transversal harmonics correspond to L = J and L — J = ±1, respectively. 
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The a abJ ^(r)) and /3£j MAr (^) are time-dependent functions and each one identifies a spin-1 mode 
from the A and Ai components, respectively. The expressions in the first equality, eqs. fll3|) and 
(|14D, represent the general expansion. In the U(l) model, we choose the homogeneous modes of 
the spin-1 field spatial components to be identically zero. In this way, the homogenous modes 
will produce field strenght configurations compatible with FRW geometries, i.e. a diagonal 
energy-momentum tensor, vanishing for the U(l) case. In addition, the r.h.s. of the same 
expressions correspond to a decomposition of the expansion in homogeneous (first term) and 
inhomogeneous modes for the non-Abelian case. There, we use a S'0(4)-symmetric Ansatz for 
the homogeneous modes of the vector field which is compatible with the FRW geometry fl6 |. 
For the case of a SO (3) global symmetry, we have, for the homogeneous modes 



Mv) = ~ 



= 0, (15) 
CMc%c, (16) 



( 2a / X 

1 + \h-xiv) 



3tt' 



where xiv) is time-dependent scalar function. The idea underlying this Ansatz for the non- 
Abelian spin-1 field consists in defining an homorphism from the isotropy group SO (3) to the 
gauge group. This homomorphism defines the internal transformation which, for the symmetric 
fields, compensates the action of a given 50(3) space rotation. Hence, the above form for 
the gauge field, where the A Q component is taken to be identically zero. By imposing the 
above mentioned symmetry conditions, we obtain a one-dimensional mechanical-type model 
depending only on time [ 1(J . The resulting one-dimensional model has some residual symmetries 



from the ones of the full four-dimensional theory. In particular, the invariance under general 
coordinate transformations in four dimensions leads to an invariance under arbitrary time- 
reparametrizations. However, due to our choice of SO (4)-symmetry conditions for the spin-1 
field, none of the local internal symmetries survive as all the available internal transformations 
are required to cancel out the action of a given SO (3) space rotation. 

Let us now turn to our model with a massive vector field. From actions (0) and (B[) one 



can work out the effective Hamiltonian density. Upon substitution of the expansions (|4])-(14D 
and after integrating over S 3 , the canonical conjugate momenta of the dynamical variables are 



found to be 



Tt„ = — = , 7T V = — = . 17 

da N dj( N 



dC eS pabNM 



3 dpfj NM 2Nna 

d£cS J f NM X /167T»J(J+1) 



(18) 



/ 1 , - — a abJ %(-lY J \ v \ (19) 

JJ d(3fj NM 2Nna 2%a N N V ; V 2J+1 ' V ; 

where C eS denotes the effective Lagrangian density associated with fl2|), which we omit here, 
and the dots represent derivatives with respect to the conformal time. To second order in 
the coefficients of the expansions and all orders in a, one finds that most of the gravitational 
degrees of freedom are gauge type and, as such, the wave function cannot depend on them; 
furthermore, we shall consider the gravitons in the ground state. For the case of SO{3) global 
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symmetry, dropping the primes, it follows that the effective Hamiltonian density reads (where 
for the Abelian casef^, one disregards the last four terms): 



n 



<-ir 



J,L 



+ [ a7T 7r ^ A L / 7r /3x? A y + Plj Pnm (L + j + iy 

\J-L\=1 

r-2 



f \, ,mj /16tt 2 J(J+1)\ 3tt 



1 



a{t) 



7T a ah J J TC ftabJJ 



+47ra 2 m 2 




abNM iqabLJ 
NM 



oabNM a 
PhJ P. 



(20) 



where, from now on, for brevity, we shall colectively denote (J,L,M,N) by J, being implicit the 
difference between longitudinal and transversal modes as well as the sum over contracted 50(3) 
group indexes. In the Appendix, neglecting the expansions (jSD-( |IO|) , we present the complete 
effective Hamiltonian for the perturbed non-Abelian model. 

The Hamiltonian constraint, 7i eS = 0, gives origin to the Wheeler-DeWitt equation after 



promoting the canonical conjugate momenta (17), (19) into operators 



7T„ 



.d_ 

da 



.d_ 

'dx 







7T aabNM 
^LJ(JJ) 



' P)fQ abLJ ( JJ ) 

°Pnm 



7T„ 



,d_ 

da 



da 



(21) 



where in (|21|) the last substitution parametrizes the operator order ambiguity with P being 
a real constant. The Wheeler-DeWitt equation is obtained imposing that the Hamiltonian 
operator annihilates the wave function $f[a, x, Plj* 1 ^, /?jj VM ]- 

Before we proceed in discussing the way the decoherence process takes place within the 
framework of our model let us comment on some features of the Hamiltonian density (^). The 
first line in (P0| ) correponds to the contribution from gravitational degrees of freedom in their 
ground state. The third line is associated with the transversal modes (J — L = ±1) of the 
spin-1 field. The fourth line corresponds, on its hand, to the contribution of the longitudinal 
modes (J = L). Notice that the presence of the mass term in the model entangles all modes 
of the spin-1 fields (longitudinal and transversal) with the metric as exhibited in the second 
line in (^). The mass of the vector field is also present in the longitudinal kinetic part, where 
there is also a term proportional to a(t) -1 . In the fifth line, in (^0]), one has the kinetic 
piece of the homogeneous part of the spin-1 field as well as the quartic potential, typical from 



3 For other presentations of an Hamiltonian formulation of systems involving the Proca field, see e.g. Refs. 

MM. 
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the dimensional compactification procedure for treating the cosmological problem of coupling 
gravity with Yang-Mills and vector fields [Q'lIIlHlIil- The coupling between the homogeneous, 
and all the inhomogeneous modes of the vector field are shown in the sixth line. Finally in 
the seventh line in (20) one has the coupling of the gravitational ground state mode to the 
homogeneous part of the mass term of the spin-1 field. Thus, our system can be regarded as 
a set of coupled harmonic oscillators (gravitational and spin-1 field) where the kinetic term of 
the longitudinal modes has a coefficient proportional to a(t) . 

In what follows, we shall compare our models with the ones discussed in the literature when 
studying correlations between coordinates and momenta and the decoherence process. 



Decoherence and Back- Reaction Processes in the Pres- 
ence of Massive Abelian Vector Fields 



In this Section, we discuss the process of decoherence in a closed FRW model with a massive 
vector field and global U(l) symmetry and its relation to the conditions which support correla- 
tions between coordinates and canonical momenta. The validity of the semi-classical Einstein 
equations, i.e. the the so-called back-reaction problem, is also discussed. As explained in the 
previous sections, we take the homogeneous modes of the spatial components of the vector 
field to be zero, whereas the environment is identified with the inhomogeneous modes, i.e. the 
/9j-functions (see eq. fl2"2p below). Therefore, our minisuperspace will be one-dimensional and 
the FRW scale factor will correspond to the physical observable about which predictions can 
be made. 

We rewrite the Wheeler-DeWitt equation (pO|),(pID for the U(l) 



H, ty[a, A A 



d 2 



2M| da 2 
j 



9n 



d 2 



d(/3 



J\2 



- Q 2 j(a) (3 j (3 j 



*[a, A^] = 0, 



(22) 



where 



an 



an 



1) 4J 48tt 3 J( J + l)/((2 J + l)4m 2 ) [1 + (l/a(t))] 



if 
if 



and 



a 



[(4m 2 a 2 /3vr) + (L + J + l) 2 ]/an 

(4m 2 a 2 /37r)+(2J+l) 2 +l 



if 
if 



J 
J 



L 
L 



J-L 
J = L 

±1 



±1 



(23) 



(24) 



I OTr+[(-l)«487r3j(J+l)/((2J+l)4m 2 )] [l+(l/a(t))] 

following the notation of Ref. M and setting the ordering ambiguity factor to P 



0. Our 

massive U(l) model and the k = +1 FRW minisuperspace model with a massive conformally 
coupled scalar field discussed in Ref. share some similar features. Indeed, up to different 
constant coefficients and the 1 + 1/ a factor, the Hamiltonian ( p0|) is equivalent to the one of a 
massive conformally coupled scalar field in a closed FRW background. As far as the 1 + 1/ a 
factor is concerned, if one considers expanding solutions, then the condition 1 + l/a — > 1 
will be rapidly satisfied and, hence, we can apply the framework used in and draw similar 
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conclusions. However, the case of contracting solutions as well as the interference between the 
two WKB branches have to be addressed differently. 



A solution of (|22)-(pl[) which corresponds to a "classical" behaviour of the a-variable on 



some region of minisuperspace will have an oscillatory WKB form (|T|) as 

tt (ft) [a, A,] = e m r s ^C {n) (a)^ {n) (a, A,) . (25) 

After expanding the functions in ( p5[ ) in powers of Mp and using (p2"|)-(p4|) one finds that the 
lowest order action, So, satisfies the Hamilton- Jacob i equation, 

- l -S% + V(a) = , (26) 



where V(a) = —a 2 + |p a 4 and the prime denotes derivative with respect to a. From (22)-(24) 



we see that the different modes do not interact among themselves. Thus, the wave function 
■0(n)(a, Afj) can be factorized as 

V>(n)(a,A M ) =^(n)(a,{/5j}) = Y[^(n)j(r],(3j) . (27) 

j 



(28) 



Defining the WKB time as 

d_ _ d^d_ 
drj da da ' 

one obtains the Schrodinger equation satisfied by each wave function ij)( n )j(r}, 

\fM (-Q§jy 2 + ^(/? J ) 2 ) An)AV, Pj) = ifywAv, h) ■ (29) 

We stress that ip( n )j(V-> Pj) * s actually ipr n )j[a(r]), /3j], dependent on the "classical" physical 
observable and f3j as well. From (p9|) , we can calculate, say, ipf n )j(T), (3j) in r( given the value 

of ip(n)j(v,Pj) at V" < i- 

In order to make predictions concerning the behaviour of a, one uses a coarse-grained 
description of the system working out the reduced density matrix associated with the WKB 



wave function of the form (Eq) to obtain 25 



PR = E ^^^^(^^(oOC^Caa^Caa, a,) , (30) 

n,n' 

where 

I n ,n'(a 2 ,a 1 ) = J ip* n/) {a 2 , A fl )^p {n) (a 1 , A^)d[A^] = Ilj J ^* J(nl) (a 2 , pj)^j( n )(ai, pj)[dpj] . (31) 

The subindex (n) labels the WKB branches. The term T n ^ n i(a 2l a\) describes the influence of 
the environment on the system. Notice that in (|3~TD all modes must be included. 

The analysis of correlations between minisuperspace coordinate and momenta is, in quantum 



cosmology, usually discussed using the Wigner function criterion 0, [30], [51]]: A strong sharp peak 
is likely to be located close to a classical trajectory defined by the Hamiltonian-Jacobi equation 
plus quantum-corrections. However, the Wigner function associated with the reduced density 
matrix ( |3~0|) does not have a single sharp peak even for a pure WKB function as (|25"D or a linear 
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combination of them (cf. Refs. [23], [32], |33], [3 



Nevertheless, this problem can be overcome 
through the environment interaction with the "observed" system |[24|| . As explained in the 
Introduction, such interaction is at the origin of the loss of quantum-coherence or decoherence 
between different classical trajectories, i.e. WKB branches. More precisely, correlations between 
coordinates and momenta must be analysed within each classical branch (n = n'). This can be 
done by looking at the reduced density matrix associated with it or the corresponding Wigner 
functional: 



dA[SU ai )SUa 2 )rh- 2 ^ A e lA ^^ (32) 



where A = (a\ — a^jl. A correlation among variables will correspond to a strong peak about 
a classical trajectory in the phase space. Thus, there exists an important relation between 
correlation and decoherence as one needs the latter, i.e. fairly small off-diagonal terms in (|30| ) 
such that quantum interference between alternative histories is negligible (X n n / oc 5 n ,n')> m 
order to obtain the former. Hence, the decoherence process is rather crucial as it is only when 
the decoherence between different WKB branches is sucessful that correlations may be properly 
predicted. 

Besides the decoherence between different WKB branches, the environment interaction also 
affects the correlations within a classically decohered branch; this is explicit in the functional 
Zn,n( a 2,ai) in eq. (^). As pointed out by Zurek the environment degrees of freedom 
continuously measure the physical observables and this interaction not only suppresses the 
off-diagonal 



(n ^ n') terms in (|30"D and ([H]), but also induces a "localizing" effect on the 
classical variables within each WKB branch. This corresponds to the back-reaction from the 
environment on the semi-classical evolution of the system. In particular, I n , n (ci2, cti) will be 
damped for |a 2 — ai| ^> 1 and the reduced density matrix associated with (^) will be diagonal 
with respect to a. The sharpness and position of the peak will be determined by the behaviour of 
1n,n{^2, 0j\). Furthermore, it has been shown in Refs. 0, ^| that the localization effect inside a 
classical branch is much more efficient than the decoherence between diferent WKB branches. 
It has been also remarked in Refs. [0, H that, if the conditions for achieving an effective 
localization (and diagonalization) of (|32| ) are met, then the interference between the different 
WKB branches is also highly supressed. Actually, the functional T n n /(a,a) has been usually 
identified as a measure of the decoherence between two different WKB histories, characterized 
by the parameters (n) and (n'); an heuristic argument in support of that view was presented 
inRef. §. 

Before proceeding, we point out that using a Gaussian Ansatz for the environment state as 

Ma,Pj) = Dj(t)e^- 



-Bj(t)pj 



(33) 



for each ipj(t], /3j), where Dj,^j are real, Bj(t) = B r j + iBij, where B r j and Bij are also real 
and B r j > 0, and imposing the normalization condition 

(34) 



the general form of X nn i(a, a') for any mode is given by 



Z(nn>)j(a, a') = exp % (D( n ,)j(a') - D {n)J {a) 



^B(nyj(a')B(n) r j{a) 

B{n')j{a') + B (n)J (a] 



1/4 



(35) 
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The analysis of the decoherence between different WKB histories and correlations via the 
Wigner function (with Z( n>n )j(a, a')) requires the functions B and D to be found explicitly. The 
requirements for sucessfull diagonalization and "localization" were generally established and 
discussed in Refs. ]7|, §. From the assumption that \a\ — <C 1 and the Gaussian Ansatz 
for ipj(a,(3j), the above mentioned conditions read: 



2M\ 



V{a) + 




\3l 



■r . \B l l\"Cf 

,7 ^-°rj 



(36) 



(37) 



(38) 



where a = (a\ + 0,2)/ 2 (the sums in these expressions and previously related ones have an 
implicit factor arising from the degeneracy of the J,L,M,N mode). Expressions (|3"(3|)-(|3"8"|) 
are usually referred to as adiabaticity, strong decoherence and strong correlation conditions, 
repectively. Notice that the results (|36|)-(|3"8|) arise directly from (|31"|),(|32|) and the Gaussian 
Ansatz for the state of the environment (1331). Henceforth, the validity of these conditions has 



to be analysed using the quantities and parameters relevant to our particular models. The 
possibility that the massive spin-1 field models give rise to new conditions for the process of 
correlation and decoherence has to be properly considered. The adiabaticity condition warrants 
the validity of the zero-th order WKB evolution as its violation implies that the semiclassical 
Einstein equations are not valid due to contributions of high-order in the phase of Z n)Tl (a2, a\). 
On the other hand, fl3"T|) reflects the fact that the peak in the Wigner function (shifted away 
from the expected classical trajectory by interaction with the environment) is sharp as far as the 
center of the peak is large when compared to the spread. Finally, expression ( |3"5| ) translates the 
condition of strong decoherence corresponding effectivelly to the requeriment of diagonalization 
of the reduced density matrix associated with (|3^). It is important to mention that usually a 
compromise between decoherence and correlation is needed since if the later is too strong, then 



the peak in the Wigner function is actually broaden [23 



Let us now address the issues of decoherence, correlations and back-reaction in our model 
with massive Abelian vector fields. Firstly, we shall assume that the decoherence between the 
two different WKB histories has occured sucessfully and consider the correlation and localization 
effects within a classical branch. Afterwards, we shall comment on the decoherence between 
different WKB branches. The inclusion of transversal as well as longitudinal modes in eqs. 
(|30D , (|31"|) , (|36|)-(p8|) give rise to some difficulties for our U{1) and SO(3) models as far as 
the longitudinal modes are concerned. This will be discussed in the following. In particular, 
notice that eq. (30)-(32) and then (36)-(38) involve considering all modes, longitudinal and 
transversal. The same will apply to other equations as will be pointed out explicitely. 

We start by considering the correlation and localization effects within each classical branch 
for the case of tranversal modes^ (J 7^ L). One can easily verify that up to a re-scaling of 
the f3 j£,-modes by a factor of mr (and conversly for the corresponding canonical momenta), 

4 The computations corresponding to each of the two WKB branches (expanding and contracting) are the 
same in the two cases (n = ±1, say). 
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this part of the model is equivalent to the one with a massive conformally coupled scalar field 
(cf. Refs. |7|, p6fl). Substituting (|3^) into the Schrodinger equation (p9|), we get the following 
equations 



7j 
Bj 



-fj(a)B rJ , 
ifM[-2{B*-tfj/4)) 



(39) 
(40) 



for Dj = 7r 1 / 4 (25 r j) 1 / 4 . With the above mentioned re-scaling, eq. ( [40]) can be linearized via 
the choice Bj = —iipj/2ipj, to yield 



(pj + Q 2 j!fj = . 



(41) 



The initial state of the environment is associated with a particular choice of initial conditions 
when solving the preceding equation. In our present case, the Hamiltonian eq. ( P2"| ) corresponds 
to a set of harmonic oscillators with a variable time-dependent mass. A convenient vacuum 
state can be defined assuming there exists an adiabatic zone such that the classical evolution 
emerging from the Hamilton- Jacobi equation satisfies the condition a/a —>■ for large a. This 
requires that V(a) be quadratic in a, meaning that models with non- vanishing cosmological 
constant do not satisfy this adiabaticity condition as, in this case, V(a) ~ 0(a 4 ) for large a. 
Notice that the Hamilton- Jacobi equation for a generic V(a) has real solutions for a > Oq only 
if ao is a single zero of V(a) for V(a) > 0. Assuming a vanishing cosmological constant, one 
can identify a vacuum for the adiabatic out regime (a ^> 1), being the out-modes of the form 



v>3 



out 



(20/)- 1/2 exp 



i / Qj(r/)drf 



(42) 



which diagonalize asymptotically the Hamiltonian for large values of a. For small values of a, a 
preferred initial state (an in vacuum state) may in some cases also be defined as the one which 
diagonalizes the Hamiltonian for a > a (in our model for A = and for other commonly used, 
da/dn <C 1, for small values of a). The relation between the in and out modes is given by the 
Bogolubov transformation 



(43) 



The aj,/3j are designated as Bogolubov coefficients and any particular choice for these de- 
termine different vacuum states for ip. One obtains for a general potential V(a) inducing an 
adiabatic behaviour [p5f 



&j - 1 ; Pj 



2 exp 



-^n[m 2 a V'(a )]~ 1/2 (d 2 j + m 2 a 2 ) 



(44) 



where the coefficient dj denotes the degeneracy associated to J, L, M, N. Notice that for the 
case of a quadratic V(a) one obtains $j = 0. If the condition dj ~ 1 ^> $j ~ is satisfied, which 
occurs when the evolution is indeed adiabatic, then such quantum state is usually identified as 
adiabatic vaccum and holds during all the evolution. Notice again that, in the case of a non- 
vanishing cosmological constant, the adiabaticity condition (Clj/fl 2 <C 1) cannot be satisfied. 

Since the massive spin-1 transversal modes behave effectivelly as conformally coupled mas- 
sive scalar fields, it comes as no surprise that the adiabaticity, strong decoherence and strong 
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correlation conditions (|36|)-(13 
Notice then that 



are indeed satisfied restricted to those environment modes. 



E 



\3l 



E 

j 



(45) 



Within the adiabaticy evolution requirement, we take a quadratic V(a) and hence the first term 
in p5|) will be dominant. Using eqs. (23) and (24), we see that the sum in eq. (^) implies (cf. 
Ref. 0) that the strong decoherence condition (|38| ) corresponds asymptotically to 



m 3 a 3 > 1, 



while the strong correlation condition (|37| ) reduces asymptotically to 

|V(S)| > m^a, 



(46) 



(47) 



for large a. As long as we restrict ourselves to assess the decoherence process and correlation 
analysis for the case of an expanding solution (as it is for the cases studied in most of the 
literature) the conditions (46), (47) are valid. In fact, that seems to be the right interpretation 
when explaning the quantum to classical transition of our Universe. As one can see, the mass 
of the of the Abelian vector field provides a scale at which the process of decoherence and the 
analysis of correlations take place. Indeed, for fairly small or negligible mass the process of 
decoherence is not completely achieved, which is consistent with the decoupling between the 
gravitational and the Yang-Mills field (m = 0) [12]. The adiabaticity condition implies that 
the growth of V(a) must be slowlier than that of a quartic potential. For a quadratic V(a) this 
is immediate. 

One could, for instance, consider an ad hoc potential from the start M, although it would 
remain to be verified if such a potential would satisfy conditions fl3~6"|)-(|3~8D and could be derived 
from a realistic action (with anisotropy or even higher curvature terms). On its turn, the pres- 
ence of a cosmological constant induces divergences in the decoherence factor (flop in addition 
to the ones from the back-reaction factor. If the latter were expected to correspond to the zero- 
point energy of the fields, the former may only be cured via the introduction of a fundamental 
cut-off since it cannot be removed by standard renormalization procedure J7|, ||. However, that 
seems rather unsatisfactory from the physical view point. Considering the environment com- 
posed by modes whose physical wavelength is larger than the horizon is in disagreement with 
the expectation that the environment consists of small wavelenght fluctuations. 

Let us now return to the decoherence between different WKB branches. The quantity 
to analyse is X nn /(a, a) (see above remarks and Ref. |§). We suppose once again that our 
minisuperspace model has an adiabatic out region for large values of a. Then, up to second 
order for the adiabatic limit and /5-Bogolubov coefficients, one can find, in the case of the 
transversal modes of the massive Abelian field IBl PSI, E26H: 



T nn f[a, a) 




-~ ( \(3 n j\ 2 + \(3n'j\ 2 + 2 (3 nJ (3 n ,j cos 



/ 

2 / (O n j — fl n / 



+ 



1 (^nj — ^n'j) 2 



16 



(4£ 



where a sum over the environment transversal modes is understood (the longitudinal modes 
will be treated in the next paragraph). The case of a quadratic V(a) imposes /3jx = and, 
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as mentioned previously, this means that the quantum state of the environment evolves as an 
adiabatic vacuum. In that case, for Q n = —Q n r we have 



Inn' (a, a) ~ exp 



1 (fl n j 
16Y 



(49) 



and, since <C 1, these terms are effectively destroyed for the two WKB histories. Another 
possibility is to consider that the Universe has undergone a static or quasi-static period for large 
a such that Cl ~ 0. In that case, the framework established above for the ipj modes can be used 
for the environment and a natural in vacuum state can be defined, which is not equivalent to the 
Associated particle creation takes place and, as a consequence, decoherence 



out vacuum 25 



occurs and the relation between the amount of interference that is suppressed and the number 
of particles created is given by 



l nn , (a, a) ~ exp | -4 ^ d 2 j[3j cos 2 



n 



nJ 



(50) 



Finally, one has to consider the longitudinal modes (J = L). Surely, they have to be included 
in eqs. (45) and (48)-(50) if the results (46), (47) and loss of quantum-interference among the 
expanding and contracting WKB branches are to be extended to a mode-complete massive 
Abelian vector field. As far as expanding solutions are concerned, namely the ones for which 
a becomes much larger than 1 sufficiently fast, the previous remarks for the transversal modes 
can be extended to this case, up to the "mode by mode" re-scaling of the /3j j functions and the 
corresponding canonical momenta. Therefore, when analysing correlations and diagonalization 
of the reduced density matrix within the WKB branch corresponding to an expanding solution 
of the Hamilton- Jacobi equation, (26), the results obtained in the previous paragraphs can 
be extended to encompass the longitudinal modes as well. The conclusions concerning the 
38|) , i.e. (45)-(47), are thus valid in the case of expanding closed FRW models 
with massive Abelian vector fields. The same applies to (48)-(50). However, for contracting 
solutions of the Hamilton- Jacobi equation, there is a problem due to the 1 + 1/a factor. The 
equation corresponding to fl4T| ) is the following 

/-I + fjj(a)\ -2 , 2 n /ri\ 

fjj + <Pjj + ttjjfjj = . (51) 

V fjj ) 

As it stands, we do not know of any exact or even adiabatic solution of this equation, which is 
of the type 

<p% = &jj<p'5f + $jj{<ft3f)*, (52) 

together with (0). It is therefore difficult to draw any conclusions about correlations as well 
as diagonalization (and localization) when the longitudinal modes are taken into consideration 
in the case of contracting WKB solutions. 

This problem becomes somewhat more acute when we try to analyse the decoherence be- 
tween the two WKB histories for n = ±1. Once again, we do not know how to obtain a form 
similar to ([48]) for X nn i (a, a) as that would involve adiabatic solutions of the type (|4^) and that 
seems, for the moment, difficult to obtain for the longitudinal modes in a WKB contracting 
solution. 



conditions (pq)- 
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A possible alternative to deal with the longitudinal modes and construct WKB solutions, 
may involve computing the influence functional (and the functions B) assuming another ap- 
proach. Namely, that the vector field mass may be almost negligible^ and therefore the presence 
of longitudinal modes could be treated using a perturbative scheme |$) , assuming that the mas- 
sive transversal modes are dominant as they are anyway present in the massless case. In Ref. 
[|7| , it was considered that the interaction between the system and the environment is such that 
there is decoherence between different WKB branches. This has been shown for specific models 
in Refs. [0, |25], [26[] for N = 1 and N > 1 minisuperspace cases. However, the use of massive 



Abelian vector field models raises the problem that for some situations that cannot be so easily 
achieved and demonstrated. Nevertheless, the factor 1 + 1/a indicates that some divergences 
will be present in the computation of the influence functional. One sould try to cure them 
by some type of renormalization procedure, eg. via a cut-off, although its nature seems to 
indicate that one has to go beyond quantum cosmological models arising from Einstein theory 
plus matter, possibly considering an effective model arising from higher- derivative theories of 
gravity or even string theories. 



4 Decoherence and Back- Reaction Processes in the Pres- 
ence of Massive Vector Fields with SO (3) Global Sym- 
metry 



In this section we extend the analysis of decoherence and correlations to massive vector 
fields with non- Abelian SO (3) global symmetry. The corresponding Wheeler-DeWitt equation 
(|20D , (|21|) can be rewritten as: 

1 d 2 

'V 



_ r _ M 2 — a 4 

2M 2 da 2 p \ 9vr 



d 2 a 
+ 



dx 2 37r 

J 



X' 



3tt 

m 

& 



+ A7ca 2 m 2 




X 



(ab) n(ab) 



Q^J(ab)^ ^J(ab)( a ) Pj P J 



tt[a,A M ] = 0, (53) 



where 
fj{ab) (a) 

and 



fi J(o6)(°) 



an + 



an if J — L = ±1 

■l) 4J 487r 3 J(J+l)/((2J+l)4m 2 )] [l + (l/o(t))] if J = L 

(54) 



{(4m 2 a 2 /3vr) + (L + J + l) 2 + 4 [l + Jf X ] /3n}/an if J - L = ±1 



(4m 2 a 2 /3jr)+(2J+l) 2 + l+4 



3tt 
2 

/3tt 



a 7r+ [(_l)4.7 487r 3j(j + i)/((2J + l)4 m 2)] [l+(l/a(t))] 



if J = L 



(55) 



5 This hypothesis should, however, be considered with some care in view of eqs. (46), (47) and their implica- 
tions regarding the presence of a massive vector field as providing a scale relatively to the decoherence process 
and correlation analysis. 
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and, furthermore, we have set the ambiguity factor to vanish, P = 0. Our minisuperspace is now 
two-dimensional, the scale factor and function x(v) (parametrizing the homogeneous modes of 
the non-Abelian massive vector fields) being the classical observable degrees of freedom. The 
environment corresponds, as before, to the inhomogeneous modes, i.e., the /^j^-functions. 

Actually, only a few particular multi-dimensional minisuperspace models have been consid- 
ered from the point of view of decoherence and correlations between coordinates and momenta. 
The Kantowski-Sachs model (N = 2) with a cosmological constant and massive inhomogeneous 
conformally coupled scalar field modes was studied in Ref. [2E], the Bianchi type-I (N = 3) with 
mass/ess inhomogeneous conformally coupled scalar field modes was studied in Ref. || and in 
Ref. PHI , an N = 2 model has been analysed where the two classical degrees of freedom were 
the scale factor and homogeneous mode of a minimally massless scalar field, the environment 
being associated with the inhomogeneous perturbations of another minimally massless scalar 
field. The analysis of mult i- dimensional minisuperspace models is, in particular, also relevant 
as it provides a possible relation between the notion of decoherence between WKB branches 
and the decoherence between histories in the the so-called Consistent Histories approach || in 
terms of space-time histories (see section V. of Ref. ||). However, no arguments have yet been 
put forward to relate the diagonalization of the reduced density matrix within a given WKB 
branch to the notion of decoherence of different histories. Nevertheless, as pointed out in Ref. 
[0, the functional X nn {a, x, of, x', •••) can be related somehow to the notion of decoherence 
between histories for the cases where N > 1. A particular solution of the Hamilton- Jacobi 
equation generates a N-l parameter family of trajectories, but there will be only one classical 
trajectory passing through each point of the minisuperspace generated by that solution of the 
Hamilton- Jacobi equation. In this sense, 1 nn {o,, Xt ■••'■> & > X > •■■) will strenghten the suppression 
of interference between histories belonging to a given WKB branch as it produces a more ef- 



ficient diagonalization of the reduced density matrix. However, calculations of Ref. |24j have 
shown that a sucessfull decoherence between histories associated to X nn has not been achieved 
for any of the models considered so far, possibly due to their simplicity |/J. 

As far as our case is concerned, we can see from eqs. (|53"D- (|55|) that the typical quadratic 
potential of massive conformally coupled scalar field with homogeneous modes is now replaced 
by the double-well quartic potential ^ x 2 ~ §= • The remaks made above concerning the 
1+1/ a factor still apply here. Notice, however, that one needs to consider the (ab)-SO(3) group 
indexes. It follows in particular, that eqs. (|27|) -(|3~2"f) (and subsequent ones) remain valid provided 



we include the SO (3) group indexes (eg. Hj — > Hj(ab)) and the x(v) function together with 
the scale factor. Generally speaking, all our working hypothesis, considerations and arguments 
presented in the last section can be extended to the non-Abelian case. The Hamilton- Jacobi 
equation is now 





\ 2 — 




3tt- 




~x 2 - 




+ 37 


25. 



+ A-na 2 m 2 




(56) 



where V(a) = -a 2 + |£ a 4 . 

As our minisuperspace is now two-dimensional, the Hamilton- Jacobi equation (|56| ) is ex- 
pected to have a one-parameter family of solutions, each one generating a family of classical 
trajectories in minisuperspace, each WKB branch interpreted as describing a whole family of 
classical trajectories, i.e. a set of different universes (and not a single one as for the N = 1 



16 



case). From the Hamilton- Jacobi, eq. (|56|), it follows that the trajectories in our N = 2 min- 
isuperspace are far more complicated than those of the N = 1 case (Section 3). Moreover, the 
WKB time is now defined as 



(57) 



= C AB dS 9 

di] ' dq A dq B 

with A,B = 1,2, qi = a,q2 = x an d Gab = diag(— 1,1). We shall have as many ^-afline 
parameters as different values of the (n)-parameter. Hence, different values of (n) will lead to 
different definitions of time for the Schrodinger equation (|29|). This implies that the influence 
functional in (|30D , (pT|) is actually a functional of two histories. The state ^J(n) (a, Xi {/3j(a&)}) 
can be interpreted, not as being simply a function of a point in minisuperspace, but instead 
as a function of the whole history, which corresponds to the only trajectory that belongs to 
the (n)-WKB branch and goes through that particular point (a,x)- Such description is fairly 
similar to the one of the Feynman- Vernon influence functional f2"7j| . 

Let us now address the issues of correlations and decoherence within each WKB branch 
and compute the relevant influence functional for our N = 2 minisuperspace model. We briefly 
describe the main framework and consider, as in Section 3, the transversal and longitudinal 

modes separately. We adopt the terminology of Ref. || here as well. Defining new variables 

i, 

(omiting the (ab) indexes for convenience) 



1 " - xa ± \ua and assuming the Gaussian Ansatz ( |3"3"D for each (3^ '-mode, we can write 



(n,n)J 



{Qa, Qa) 



cxp 



A R 

-e {J) y A y B 



exp 



{hi 



where e AB and e A are designated as decoherence matrix and phase vector respectively, and 

1 



AB 

e (J) 

~A 
e (J) 



h 



[(B rJ y A (B rJ y B + (B tJ y A (B u y B } 



{Dj VA 



(B u y A /A(B rJ 



J A 



(59) 
(60) 



Let us again consider the minisuperspace with a region, say for large values of a, for which 
an adiabatic solution of eq. ((41]) of the type (pE2| ) and ([|3|) can be established. The quantum 
state of the environment is determined by choosing the Bogolubov coefficients. We then obtain 



AB 

e {J) 



~A 
e {J) 



>A/Br,2 i 



JA 



2a 



4a 



i A 



(61) 
(62) 



The first term in ( |5"T| ) corresponds to the adiabatic vacuum contribution (f3j = 0) while the sec- 
ond one is related to particle creation. It is interesting to notice that for N = 2 one can always^ 
generate a diagonalization along two independent directions of minisuperspace by coupling to 
a variable mass. 



6 For N > 2 the situation is however, different because e AB is positive definite and also degenerate whenever 
N > 2 as the directions of the (B r j)' A and (Bij)' A vectors diagonalize the reduced density matrix and therefore 
any other orthogonal direction to these is an eigenvector with null eigenvalues. 



17 



Correlations between each minisuperspace coordinate a and \ and their canonical momenta 
can be analysed by examining peaks in the reduced Wigner function 



which is equivalent to 



F wl (q A ,7T qi ) = J dn qa F W ( ri )(q A ; ir qA ] 



(63) 



F wl (q A ,ir qi ) = exp 



11 



dS 



77(11 ~ M ' P dq~ l ~ ' 



~1 



(64) 



where ir qi is the momentum conjugate to q A , A — 1. The strong correlation condition translates 



as 



C AA 



(65) 

(where iT~X is a typical value of the momentum along the trajectory) and the strong decoherence 
condition is 

e AA > l/q A . (66) 

Notice that eqs. (|63"D 



must include the summation over the SO (3)-group indexes and all 



modes, using eqs. 



to 



The same holds for eqs. 



when analysing decoherence 



between different WKB branches. 

The above construction is valid for the transversal modes as those behave similarly to 
massive conformally coupled scalar fields. Let us take the case for which j3j = 0, i.e., we choose 
the quantum state of the environment modes to be the adiabatic vacuum. In this case we must 
put A = 0. The relevant quantities to analyse the decoherence and correlation will be e 



,22 



j 1 and 



ui 



/ 4m 2 
1 3tt 



Am? a? 
3tt 



+ (L + J+ l) 2 + 4 



1 + 



— Y 

3vr A 



(67) 



22 



2a 
3?r 



1 + 



— y 



Am? a? 
3?r 



+ (L + J+ 1)2 + 4 



1 + 



-x 



(68) 



The sum in J implies that, similarly to Ref. [0, e 11 and e 22 behave, for large a, proportionally 

Notice that one expects the v-field to evolve towards one of the 
for large values of a, following Ref. [16], and neglecting 



to a and a 1 , respectively, 
minima of the potential ^ 



curvature terms as a —> oo. From the Hamilton- Jacobi equation we obtain that a typical value 
along a WKB trajectory for the canonical conjugate momenta to a and x are — ff ~ a* 
and 7r, 



Such 



, , - || ~ x Z i assuming a small mass in order to neglect the interaction terms, 
approach has already been discussed in Section 3 regarding the longitudinal modes. As far as 
WKB expanding solutions are concerned, conditions (65) and (66) with (67) and (68) seem 
to indicate that, for the \- field, the strong correlation condition will be satisfied but not the 
strong decoherence one. For the scale factor, the strong correlation and strong decoherence 
conditions will be satisfied in the very sense of Ref. |]]. We shall discuss the apparent failure 
in fulfilling the strong decoherence condition for the \- field in Section 5. 

Once again, the longitudinal modes are well behaved concerning diagonalization and corre- 
lation within a suitable expanding WKB branch and the above results are equally suited here. 
However, the difficulty associated to contracting branches still remains. 
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Equations (48)-(50), including the (ab) — SO (3) indexes, are valid for the transversal modes. 
The same conclusion holds for longitudinal modes within expanding WKB branches. If we 
consider two different WKB expanding branches in our N-2 dimensional minisuperspace with 
an adiabatic vacuum, eq. (48) can be reduced to 



X nn i (a, a) 




1 (VS (n) -VS K) )(Vfi 



64 



n 



(a6)4 



(69) 



with VS l ( n )V = G AB j^-^. Hence, we find that the interference between terms with different 

(n) is exponentially supressed, depending on how different is the WKB time variation of 
along the two trajectories (cf. comment after eq. (57)). However, when there are two branches 
and one corresponds to an expanding solutions while the other to contractiong ones, the issues 
raised in Section 3 with respect to longitudinal modes equally will apply in this case as well. 
Namely, one cannot use (48)-(50) unless some consistent perturbative scheme to treat the 
longitudinal modes when m <C 1 is available and justifiable. 



5 Conclusions and Discussion 

We have discussed the quantum cosmology of a massive vector field coupled with gravity 
and we have shown that the resulting model possesses interesting properties in what concerns 
the decoherence of the scale factor of a closed FRW geometry. In the presence of a massive 
vector field with U(l) global symmetry, the scale factor is the only decohered quantity while 
for the non-Abelian case with SO (3) global symmetry, the scale factor and the homogeneous 
mode x{t) ar e the decohered variables expected to behave classically. 

As far as we consider expanding semiclassical solutions, the models we propose can be re- 
garded, to a certain extent, on the same footing as the ones where decoherence of degrees of 
freedom of the metric is achieved via tracing out higher modes of self-interacting scalar fields. 
The inhomogeneous modes of massive vector fields, which were expanded in spin-1 hyperspheri- 
cal harmonics, represent an interesting alternative to play the role of environment for the metric 
in the Abelian case and for the metric and the homogeneous mode of the non-Abelian massive 
vector field. However, in the latter case, we find that the strong decoherence condition, eq. (66), 
for the x-field (parametrizing the non-Abelian massive vector field homogeneous modes) is not 
satisfied. Unfortunately, since the literature on iV > 1 minisuperspace models is rather scarce, 
we could not contrast our results with the existing ones. We can mention nevertheless, that, 
for instance, in Ref. [8], in a diagonal Bianchi type-I model where the environment consists 
of modes of a massless conformally coupled scalar field, a problem of similar nature to ours is 
encountered. One could speculate whether the decoherence and correlation conditions would 
be satisfied using solutions that already account for the back-reaction, rather than the classical 
histories or, instead, considering higher- derivative terms as done in Ref. ||. A self-consistent 
approach along these lines could be seen as a way to implement our massive non-Abelian vector 
field model. 

We have also shown that, for the transversal modes, known techniques and the associated 
discussion on decoherence and localization within WKB branches, namely using the Gaussian 
Ansatz for the wave function and computing the Wigner functional in the adiabatic limit, still 
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holds for both cases we have considered (Abelian and non-Abelian) . The same can be said 
for the longitudinal modes at least in what concerns expanding models. Issues related with 
the problem of treating the longitudinal modes in contracting models, in particular regarding 
the decoherence between different WKB branches, have been discussed and some possible fixes 
have been sugested. 
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Appendix 



The full effective Hamiltonian density is given by: 
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where D N is a representation for the scalar harmonics 

QX» = f?(x)*WM) , 



e d i+1 (cosnx) 



d(cosx) i+1 

are Fock harmonics, Y em (9,(f)) are spherical harmonics on S 2 and ^ are 3 — j 

symbols. y s3 g denotes the square root of the determinant of the metric over the unitary 
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3-sphere and L l a represents the transformation between the left-invariant basis on S* 3 and a 
coordinate basis. We have used 
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Furthermore, we have made use of the relation provided by the equation of motion: 



i.e., 
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(d*F l0 ^ +m 2 A ^) + ([A\F t0 }) (bc) = 0, 



where ( )^ means the (6c)-component projection. Using the expansion of A = A^uj^T^ in 
the above equation, multiplying by D J M N and integrating over S 3 we get after integration by 
parts (cf. ref. [36]): 
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From the last term of the last equation we get integrals of the type Si, S2, S3, S4, the first 
two terms are valid only for the Abelian case (where the (a6)-SO(3) group indices have been 
obviously disregarded) . 
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